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Abstract 

We describe an algorithm for converting one bipartite quantum state into another using only local operations 
and classical communication, which is much simpler than the original algorithm given by Nielsen [Phys. Rev. Lett. 
83, 436 (1999)]. Our algorithm uses only a single measurement by one of the parties, followed by local unitary 
operations which are permutations in the local Schmidt bases. 

1 Introduction 

Consider the case where two parties, Alice and Bob, share an entangled state of two JV-level particles. M. A. Nielsen 
introduced in [|N] an algorithm for converting one such pure bipartite state into another using only local operations 
and classical communication. He gave an explicit condition for when such a conversion is possible: We write the 
first state in Schmidt form 
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\ci>) = ^2Vrf\i)®\i), (i) 

i = l 

where the a 1 are non- negative real numbers satisfying a 1 = 1, and similarly for the target state 

N 

w^E^i^B- (2) 

Then \(f>) may be converted to \ip) iff the vector (ft) majorizes (ct), i.e., 

i=i i=i 

with equality for k = N; here ^ ft are the ft arranged in descending order, and similarly for the ^«\ This is again 



equivalent (theorem II. 1.10 of |Bh|) to the existence of a doubly stochastic matrix D such that a = Df3 (a doubly 
stochastic matrix has non-negative real entries, and all its rows and columns sum to one). 

Nielsen's algorithm uses several rounds of individual measurements and classical communication. Alt hou gh it 



is known that such a sequence of operations can be replaced by one involving only a single measurement [LP], the 
proof of this is non-constructive and not easily applied to Nielsen's algorithm. Our simpler algorithm may be useful 
for practical applications and for the analysis of local pure-state conversion in quantum cryptography ]b|, |jS[ . A 
similar result using a different method has been obtained earlier by Hardy Q. 



2 Example 

In this section we illustrate how the algorithm works by an example. We consider the case ft 1 = (ft) T = 
(3/5, 3/10, 1/10) and oft = (ct) T — (2/5, 1/4, 7/20) (for typographic reasons, we write the transpose of the column 



1 



vectors); note that a is not sorted: as we s hall see, this doesn't matter. We check that /3 y a. Using the algorithm 



from (the proof of) theorem II. 1.10 in [Bh , we find a doubly stochastic matrix that maps j3 to a: 



( 1/3 2/3 \ / 3/5 \ / 2/5 \ 

1/6 1/3 1/2 3/10 = 1/4 (4) 
\l/2 1/2 / \1/10 / \7/20 / 

From this matrix we now derive the set of unitary transformations, which turn out to be permutations. We start 
by finding a set of non-zero entries with no two entries in the same row or in the same column, i.e., corresponding 

(° 1 °\ 

to some permutation matrix. We first choose positions 1 I , corresponding to the permutation (12). The 

\CJ l) 

smallest entry in the doubly stochastic matrix in these positions is 1/6, so we subtract 1/6 times the permutation 
matrix. Then we get 

/1/3 1/2 \ 

1/3 1/2 (5) 
\ 1/2 1/3/ 

Again we look for a set of non-zero entries with no two in the same row or column; let's pick the identity matrix 
this time. The smallest entry on the main diagonal of Q is 1/3, so we subtract 1/3 times the identity matrix, and 
finally we are left with 1/2 times a permutation matrix corresponding to the permutation (132). 
Using the above decomposition, we can now rewrite (M) as 




1/6 \ 


/ 2/5 


1/3 = 


1/4 


1/2/ 


\7/20 



/3/10 3/5 

3/5 3/10 1/10 | | 1/3 | = | 1/4 | (6) 
\ 1/10 1/10 

We may imagine that the columns in the matrix are indexed by the permutations we found above, i.e., (12), (), 
and (132), respectively. The entries of column a are then f) permuted by cr" 1 ; we permute the vector by shuffling 
the components: (J3 1 , /3 2 , /3 3 ).(123) = (/3 2 , /3 3 , /3 1 ) (this corresponds to the last row because (132)' 1 = (123)). Thus, 
in the column corresponding to a, the ith row has the entry /3 CT The LHS vector in (|^) is (p a ), where the p" 

is the weight corresponding to a that we found above when decomposing the doubly stochastic matrix. The RHS 
vector is, of course, a. 

Now we need to find the POVM. To do this, we find three diagonal matrices A a , each defined by taking column 
a and dividing the ith entry (in the ith row) with a 1 , and multiplying with p" . Thus 

A (12) = |diag(%^,f^|, i^|§) =diag(l/8, 2/5, 1/21) (7) 

A Q =Idiag(|f,^,iZ§) =diag(l/2, 2/5, 2/21) (8) 

A (123) =|diag(^,i^,^) =diag(3/8,l/5,6/7) (9) 

We observe that each A a is a positive matrix, and that the sum of the A a is the identity matrix. The A a thus 
define a POVM. All we now need to see is what happens when the corresponding operations are applied to the 
state described by a. 

The POVM is performed locally by one of the parties, say Alice. In our example, the post-measurement state 
corresponding to the outcome a is a pure state with Schmidt coefficients obtained by normalizing the vector A a a: 

• (12): (1/20, 1/10, 1/60) ~> (3/10,3/5, 1/10) 
. (): (1/5, 1/10, 1/30) ~» (3/5, 3/10, 1/10) 

• (123): (3/20, 1/20,3/10) ~» (3/10, 1/10,3/5) 

where the ~» indicates the normalization. The post-measurement states are thus given by the columns of the 
matrix in equation (^). To reach the target state, Alice permutes the bases according to the permutation <r, and 
she communicates a to Bob who then performs the same permutation on his basis. 
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3 Formal Description and Proof 



3.1 Description of the algorithm 

Find the doubly stochastic matrix 

Given vectors of Schmidt coefficients a and /3 such that 13 >- a, we first use the algorithm of |Bh|, II. 1.10, to 
find a doubly stochastic matrix D such that a = D(3. 
Decompose the doubly stochastic matrix 

The Birkhoff-von Neumann theorem (]Bh[, chapter 2) states that we can find permutations a G E C Sjv and 
positive numbers p a such that 

where P„ is the permutation matrix associated to a, i.e., it maps the fundamental vector (with a one in the ith 
place and zeros otherwise) to e a u). 

To explicitly find the decomposition of the density matrix, we may use that many of the proofs of the Birkhoff- 
von Neumann theorem are constructive. The approach we sketch here is similar to the proof in jp| . 

It will be useful to define an arrangement (or diagonal) of a n x n matrix M — (to}) as a set of entries {m^^} 
for some permutation a £ 6 n , i.e., where no two entries are taken from the same row or the same column. It 



follows from the Konig-Frobenius theorem ([Bh|, chapter 2) that any doubly stochastic matrix has an arrangement 
with all entries non-zero. 

There is a polynomial-time algorithm for finding such arrangements in a matrix: Consider the matrix D which is 
D with each non-zero entry replaced by 1. Let R = {ri , . . . , r n } be the set of row indices of D and C = {ci , . . . , c„} 
the set of column indices. Then D defines a balanced bipartite graph G with vertices R U C (disjoint union): we 
have an edge from r; to c, iff D5 = 1. By the Konig-Frobenius theorem, G has no isolated nodes. We can thus 
easily find a perfect matching in the graph; there are polynomial-time algorithms for doing this, see e.g., chapter 3, 
section 3 of Q. A perfect matching defines a unique permutation a £ & n through its edges r a ^ <-> Cj, i = 1, . . . , n. 
It is easy to see that the matching also defines an arrangement of D with all entries non-zero. 
Constructing the POVM 

Now we assume that we have the decomposition (Jic)) . Hence 

Vi: ^/T~ 1( V = a (11) 

a 

Indeed, a 1 = (Dpy = p CT P CT )/3] 1 = E.P" ■ 
Next, we define matrices 

^ CT =p CT dia gl (^— ) (12) 

cr 

where diag i (/(i)) denotes a diagonal matrix whose entry is /(«)• 

The {A a | er € E} define a POVM, since each A a is a positive matrix, and 



^(^E^ CT " 1(i) ) 



Furthermore, 



= diagi 1 

tr^diag^a 1 )) = p" (13) 



so that p" is the probability of outcome a. 
The measurement 

Now we turn to the measurement itself. When Alice performs the measurement on her side, we need only 
consider her reduced density matrix, p, which has eigenvalues a 1 . We choose the operations corresponding to 
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the POVM A a such that, if the outcome of the measurement is a, then the reduced density matrix after the 
measurement is 

Pa = — vXpVX' = —A a p. (14) 
p p 

The eigenvalues of p a are w , i = 1,...,JV, which sum to one. Alice then performs the permutation a on her 
side, and communicates a to Bob so that he can do the same: since Y^i I*) ® !*) = IZj /^ I | cr (*)) ® l cr (*))i the 

operation should map \<r(i)) to |z) for all i. After the unitary operations, the particles are in the pure state \if)) (0) 
with Schmidt coefficients f3 l . 



3.2 The converse 

We now turn to the converse: given a and f3 and a POVM {A a \ a £ E} with £ C S n such that for all a £ E, 

(A CT a).CT =p CT /3, (15) 

for some positive constants p a , we show that /3 >- a (of course, this also follows directly from Nielsen's theorem 
|n[|). We assume that each A a is given by a diagonal matrix. 
Define a matrix Y — (7]) by 

7j= E p CT ( 16 ) 

{<T|o-(j) = i} 

Then we have 

v i ; E^ = E Y p ° 

= £y = i; 

ctSE 

Vi: Y^i = E E 

i i {o-|ct(j)=0 

ctSE 

so r is doubly stochastic. Now equation ([[5]) implies A a a — p CT /3.(o- _1 ), which again implies that A a must be of 
the form (|l|). Finally, 

Vi: (r/?) 1 = ^ 7 j^ 

i 

= E E ^" 1(4) 

3 {o-|o-(j)=*} 
ct£E 

= a*. (17) 

Thus r/3 = a, which means that /3 >- a. 
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